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Abstract 


Double-Gate (DG) MOSFETs are becoming the most promising device structures for 
scaling in nanometer range. In view of anticipated role of DG MOSFET in digital, 
analog and mixed signal device applications, modeling and simulation of the device 
behaviour has assumed considerable importance. In this work Sub-threshold behavior, 
which determines the switching characteristics and the short-channel immunity of the 
device, has been studied. Physics based analytical sub-threshold model for un-doped 
DG MOSFET is proposed. The model results were compared with numerical 
simulation results of 2-D device simulator (ATLAS). A good agreement was found 
between the proposed model and the numerical simulation results. The strength of the 
model is that it is physics based and simple and yet capable of predicting the device 
behaviour reasonably accurately. 
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Introduction 


Chapter 1 


1.1 Overview 

Bulk-MOSFET and SOI-MOSFET have been in the area of VLSI for past many 
years. Devices have been scaling down, as the process technologies are improving, for 
faster speeds. But as one is scaling down into the nanometer range the gate tunneling 
currents and short-channel effects are coming into picture, which affect the proper 
functioning of the devices. This has lead to increased research towards coming up 
with newer devices with improved performance without adding to processing 
complexity. One of the intuitive methods of controlling the short-channel effects as 
one goes down in length is to somehow increase the gate control over the channel 
length. For this, multi-gate MOSFETs like double-gate (DG) MOSFETs, n-gate or 
Fin-FET, and gate all around (GAA) MOSFET [1],[2],[3],[4] are being studied. Of 
these newer structures, double-gate (DG) MOSFET is becoming major area of 
research, as it has shown the ability to scale down to very small lengths with 
acceptable short-chaimel effects without adding significantly to processing 
complexity. It has been anticipated that it will not only be used in digital domain but 
will also be extensively used in analog and mixed signal applications. 

All around the world effort is being put to come up with a compact model for 
DG MOSFETs so that its characteristics can be analyzed by simulation and also to 
determine the maximum scaling limits to which it can be scaled down. Many models 
have been proposed till now. Some are based on the modeling of inversion charge 
[5],[6]. In the paper by Jin He et. al [5], they have started with one dimensional 
Poisson’s equation and instead of solving it for potential, they have solved it for 
inversion charge and then used that expression of inversion charge to evaluate 
potential and current in the device. But models based on inversion charge have an 
inherent disadvantage that the moderate inversion regime of device operation is not 
that effectively modeled, as is modeled by surface potential based model. For DG 
MOSFETs several surface potential based models have been proposed [7], [8], [9]. 
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One such model was proposed by Yuan Taur [9] in which he has came up with a 
model for drain current of DG MOSFETs with undoped body. But this model has a 
problem that it is not analytical model and needs to be solved iteratively. In another 
model by Malobabic et. al [10], they have simplified the model for drain current to 
come up with an analytical model. But all these models deal with long channel 
devices. Thus a model is needed for short channel devices. 

Since these devices have shown good immunity to short-channel effects and 
due to their anticipated use in analog and mixed signal applications, sub-threshold 
region becomes an important area to be given special attention. Till now many works 
have been published on DG MOSFET. But as far as sub-threshold region specifically 
is concerned, only one model for sub-threshold is proposed by X. Liang et. al [1 1] .In 
this work they have solved the 2-D Laplace’s equation in terms of series solutions and 
then they came up with an expression for the characteristic length for scaling of 
devices and also an expression for potential variation in two dimensions. The model is 
able to predict the device behaviour in sub-threshold region quiet accurately along 
with the threshold voltage roll-off, drain induced barrier lowering (DIBL) effect and 
sub-threshold slope. But the inherent problem with the model is that it is not 
analytically solvable and even after making simplifying assumptions numerical 
iterations can not be avoided. In the present work, attempt is made to come up with an 
analytical model for sub-threshold region with reasonable accuracy and much simpler 
formulation, ultimately resulting in faster calculations. 

In earlier works on sub-micrometer MOSFETs by Zhi-Hong Liu et. al [12], a 
model for surface potential in sub-threshold region has been formulated based on the 
application of Gauss’s law to the silicon body. They have taken the effect of lateral 
field also while applying Gauss’s law. Since the formulation was done for sub- 
threshold region, the inversion charge was considered negligible compared to ionized 
dopants and was neglected. They have also given the expression for characteristic 
length for scaling with only one fitting parameter in the model. The model has shown 
to predict the results reasonably accurately for lengths down to 100 nm. 

A model based on similar formulation obtained by application of Gauss’s law 
to the fully depleted SOI MOSFETs has been proposed by Srinivasa R. Banna et. al 
[13] . In this model they have applied Gauss’s law to the box of SOI MOSFET, again 
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taking into account the flux due to lateral field and neglecting the mobile charge 
carriers compared to ionized dopants. Similar expressions for simple bulk-MOSFETs 
were obtained. Again the model was able to predict the surface potential variation in 
the device very accurately. 

As already mentioned, in the case of DG MOSFETs only one model for sub- 
threshold region has been reported. The inherent difficulty in DG MOSFETs is the 
variation of potential in the direction perpendicular to channel, i.e. from one gate to 
another gate. Thus this results in need for non-analytical solution for the potential at 
the center of the silicon body film. This renders the solution to be solved by using 
numerical / iterative methods. This is true for Liang and Taur [11] model which is 
solvable only by numerical iterations. Thus there is a need for making simplifying 
assumptions to come up with an analytical model for potential variation. 

The basic device structure that is used for the formulation of models is as shown 
below. 



Figure 1.1. Basic Device Structure 

DG MOSFET can operate in two modes: Symmetric and Asymmetric Mode. In 
symmetric mode, the line of symmetry of device goes through the center of the silicon 
body or one can say that everything is same for front gate and back gate namely. 
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oxide thickness, oxide material and gate work -function. Same gate voltage is applied 
to both the front and the back gate. The device can be operated in asymmetric model 
by driving front and back gate independently or by having different oxide thickness or 
by having different gate work functions or by having combination of any of these 
methods. Thus from the above discussion it becomes clear that to come up with a 
model for such a device, both modes need to be considered separately. 

1.2 Organization of thesis 

This work has been divided into two parts based on operation mode discussed above. 
Part I will look into the formulation of the model for symmetric DG MOSFET with 
chaper 2 dealing with formulation of surface potential expression along the channel 
and chapter 3 dealing with the formulation of the sub-threshold current model, sub- 
threshold slope formulation and the comparisons of the results obtained using model 
with 2-D device simulator ATLAS [14]. 

Part II consists of chapter 4 and 5 with chapter 4 dealing with the formulation 
of surface potential expressions and chapter 5 dealing with formulation of current 
expression in sub-threshold region for asyimnetric DG MOSFET in similar way as 
done for symmetric case. 

In the end in chapter 6 a comparison is done between syimnetric and 
asymmetric mode and the models for symmetric and asymmetric modes are combined 
to give a generalized model for DG MOSFETs in sub-threshold region. 
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Chapter 2 

Modeling of Surface-Potential along the channel in the Sub-Threshold 
Region using Gauss’s Law 


2.1 Gauss’s Law 


Gauss Law is one of the most widely used laws in electrostatics. It is not only 
used in semiconductor devices but is used anywhere there is an electric flux and 
accumulation of charge. For geometries of sufficient symmetry, it simplifies the 
calculation of the electric field. In mathematical form it can be represented as. 



-> 

A = 


Q 


Ej-Sq 


where E is the electric field d A is the normal to the surface and Q is the net charge 
contained within the surface. 

In the case of semiconductors there can be two kind of charges - one due to the 
ionized dopants and other due to mobile charge carriers i.e. electrons and holes. 
Gauss’s Law will be applied within the silicon body of the DG MOSFET, which for 
present case is undoped. The device structure that will be considered for analysis is 
shown in Figure 2.1. Applying the Gauss’s Law to the silicon body of the DG 
MOSFET 


G/2 


'SI 




d(p(x,y) 
dy 


dx -2 Sox \ Eqx (y' )dy'= Qi tsi y 


( 2 . 1 ) 


where the first term on the L.H.S represents the net electric flux entering the silicon 
body from the lateral direction and the second term represents the net electric flux 
entering the body firom the top and bottom. The R.H.S represents the net charge 
within the body. 




Figure 2.1. Symmetric Double-Gate (DG) MOSFET structure 

Sub-threshold region will be considered for the present case. Hence the mobile charge 
carriers will be neglected. Thus Eq. (2.1) becomes 

( 2 . 2 ) 


-t./ / 2 


2.2 Assumption for flux entering from lateral direction 

In Eq. (2.2), as mentioned earlier, the first term on L.H.S represents the net 
electric flux entering the silicon body from the lateral direction i.e. source and drain 
side. It has been shown in the earlier works of Suzuki et. al [15] and supported by 
Yuan Taur [16] that in the sub-threshold region, the difference between the surface 
potential potential at the centre of the silicon body is very small or 

in other words it can be said that the variation of potential in the x-direction ( or 
perpendicular) direction is negligible. To support this argument figure 2.2 and figure 

2.3 show the variation of surface potential and centre potential with respect to gate 

voltage, F(j. 






Figure 2.2, Surface potential and centre potential comparison with silicon body thickness, tsi 
varying from 5nm to 20nm in steps of 5nm, oxide thickness, to* = 1.5nm and VFB,r~ VFB,b = 0.0 V 



Figure 2.3. Surface potential and centre potential comparison with oxide thickness, to* varying as 
l.Onm, 1.2nm, l.Snm and 2,0nm, silicon body thickness, t^i = 5nm and VFB,f ~ ” 0.0 V 
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Figure 2.2 and figure 2.3 are obtained by solving numerically the following equation 


^FB~ ^centre' 


r //iM 

■— ln[cos(^)J-^i-^ ^exp 


^ Wee., re ^ 

2kT 


tan (6') 


(2.3) 


where 0= 


q\ 

' 2e,,kT 


exp 


' Ween, re \ si 

y 2kT J 2 


The derivation of Eq. (2.3) is given in Appendix A. Here the surface potential in terms 
of centre potential is given as 

^s = (Peemre ^ [^ 0 ® (^)] ( 2 - 4 ) 

Thus from figure 2.2 and figure 2.3 it is clear that in the sub-threshold region the 
surface potential and the centre potential are almost the same. The results have been 
verified by 2-D device simulator [14]. Hence <p{x,y) can assumed constant in the x- 
direction. Based on this the following assumption is made 

'"P (2.5) 

i dy T] dy 

where is the fitting parameter [12], [13] which depends upon the oxide thickness, 
silicon body thickness and channel length. For short channel devices potential is not 
constant along the perpendicular direction due to strong effect of lateral fields. Still 
the difference between the surface potential and centre potential is not significant and 
hence the parameter t] has been introduced. Thus from Eq. (2.5) it is clear that the 2- 
dimensional (2-D) problem effectively reduces to a 1 -dimensional (1-D) problem. 

Eq. (2.2) and Eq. (2.5) give 


^sihi depsiy) 

V dy 




y 

]Eoxiy)dy'=^ 


0 


Now (y ) is the field coming from the gates and is given by 

p (r\ \EG-EFB-<Psiy)] 

^ox 


( 2 . 6 ) 


(2.7) 


^ Refer section 3.5.1 



Putting Eq. (2.7) into Eq. (2.6) and differentiating with respect to y 
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F(Psiy) 

dy^ 


- \ps (>") ~ (^G ~ ^FB )] 
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( 2 . 8 ) 


r— — 

where X = - Characteristic Length 

2.3 Formulation of boundary Conditions 

Now to formulate the boundary conditions an undoped body DG MOSFET is 
assumed with n+ source and drain and since the silicon body is undoped, free charge 
carriers are assumed to be negligible. 



Figure 2.4. Boundary conditions for Short-Channel Symmetric DG MOSFET 

based on Gauss law. 

The boundary conditions, as shown in figure 2.4, are 

(Ps{0) = (p{xSy)= hi 
(PsiL) = (pix,L)= hi + '^DS 


(I) 

( II ) 


'a + -'^(|) = (1'x)(I) 



Now (I) comes from the fact that source/drain is heavily doped and body is undoped, 
hence treated as a step-graded junction. The body being undoped results in a potential 
drop in the body, known as the built-in potential and can be expressed as, [1 1] 

£ 

4>bi ~ Work-function difference between n’*’ source/drain and intrinsic body. 

2q 

(II) is same as (I) only the drain junction is reversed biased by an additional bias of 
Vds and body being undoped, this potential also drops across the body resulting in 
boundary condition (II). 


2.4 Analytical Model for Surface Potential in Sub-threshold 
Region 

Eq. (2.8) is a standard differential equation and can be solved easily along 
with the boundary conditions (I) and (II)^ to give the following expression for the 
surface potential along the channel in the sub-threshold region 






(2.9) 


where - Vq- V pg 

From Eq. (2.9) it is evident that for short-channel devices, surface potential is not 
constant as is also verified in [12], [13]. Figure 2.5 and figure 2.6 show the calculated 
results based upon Eq. (2.9) along with its comparison with the 2-D numerical 
simulations obtained using Atlas device simulator. 

From figure 2.5 and figure 2.6 it is clear that as the channel length is 
decreased, potential no longer remains constant along the channel. Although the 
potential is not constant near the source/drain region but that is due to the source-body 
and drain-body junction and not due to the short-channel effects. In figure 2.7 and 
figure 2.8, the variation of surface potential along the channel for larger values of 
drain voltage is shown, showing the validity of the model for higher Vds- 


^ Refer Appendix B 


in 




Figure 2.5. Comparison of Surface Potential along channel for analytical model Eq. (2.9) and 
2-0 Numerical Simulation (short lengths) for small drain bias (Vds “ 0.01 V)» to* == 1.5 nm, 
t„ » 5nm, Vg » 0.1 V, = VFB,b =* 0.0 V, L*8, 10, 20, 30 and 50 nm. 
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Figure 2,6. Comparison of Surface Potential along channel for analytical model Eq. (2.9) and 
2-D Numerical Simulation (longer lengths) for small drain bias (Vds = 0,01 V), t„ = 1.5 nm, 
t,, = 5nm, Vc = 0.1 V, Wfbj = = 0.0 V, L=70, 100, 200 and 300 nm. 
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Figure 2.7. Comparison of Surface Potential along channel for analytical model Eq. (2.9) and 
2-D Numerical Simulation (short lengths) for large drain bias (V^s === E5 V), t^x = 1.5 nm, 
t,j * 5nm, Vg = 0.1 V, VFB,r*' VFB,b 0.0 V, L*8, 10, 20, 30 and 50 nm. 
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Figure 2.8. Comparison of Surface Potential along channel for analytical model Eq. (2.9) and 
2-D Numerical Simulation (longer lengths) for large drain bias (Vns “ V), t„, = 1.5 nm, 

t., = 5nm, Vg = 0.1 V, VKB.f = Vre,, = 0.0 V, L=70, 100, 200 and 300 nm. 
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In figure 2.9 and figure 2.10 both plots of surface potential for low and high values of 
Vps ^re compared and it clearly shows that for long channel devices even for higher 
drain biases, the minimum potential along the channel remains unaffected while for 
shorter length devices, the surface potential is raised clearly portraying the effect of 
lateral field. This increase in minimum potential along the channel due to increase in 
lateral field is effectively reducing the barrier to mobile charge carriers to flow from 
source to drain and hence results in reduction in the threshold voltage of the device. 
This effect of reduction of threshold voltage with increase in drain voltage is known 
as Drain Induced Barrier Lowering (DIBL). Even at lower drain biases, the minimum 
potential along the channel keeps on increasing as the length is decreased as was 
shown earlier in figure 2.5. This is due to reduction in channel length. Thus is known 
as Short-Channel Induced Barrier Lowering (SCIBL) [18]. This effect of reduction in 
threshold voltage with decrease in channel length has been attributed to increase in 
the charge shared by the source-channel and drain-channel depletion regions. This 
effectively reduces the amount of charge required by gate to form the channel and 
hence the threshold voltage. 



Figure 2.9. Effect of drain induced barrier lowering (DIBL) for short lengths and comparison 
with 2-D Numerical Simulation t„ = 1.2 nm, t„ = Snm, Vc = 0.1 V, VFB,r= VFB,b = 0.0 V, L=8, 10, 
20, 30 and 50 nm and Yds “ 0.01 and 1.5 V 
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Figure 2J0. Effect of drain induced barrier lowering (DIBL) for short lengths and comparison 
with 2-D Numerical Simulation tox = 1.2 nm, t,i = 5nm, Vg = 0.1 V, VpB.r = VFB,b = 0.0 V, L=70, 
100, 200 and 300 nm and Vos = 0.01 and 1.5 V 

2.5 Effect of Oxide Thickness on surface potential 

Using Eq. (2.9) surface potential is plotted along the channel and compared 
with the 2-D numerical simulation showing the variation in surface potential due to 
the variation in oxide thickness. Figure 2.11 gives a clear picture that keeping 
everything else constant and reducing the oxide thickness results in reduction of 
minimum value of the surface potential along the channel showing that more gate 
voltage will be required to form the channel for thinner oxides and thus predicting 
increased threshold voltage with decrease in the oxide thickness. This effect is more 
pronounced in the shorter length devices compared to long channel devices. Similar 
things were observed by Wong et. al [19]. 

2.6 Effect of Silicon Body Thickness on surface potential 

Looking at figure 2.12 it can be found that as the silicon body thickness is 
decreased the minimum potential along the channel keeps on decreasing and thus 
showing increased coupling between the two channels formed in the device. This 
results in increased volume inversion in the silicon body. Thus more gate voltage is 
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required to form the channel under the gates, resulting in increase in threshold voltage 
with decrease in the silicon body thickness. 



Figure 2.11. Effect of oxide thickness variation on surface potential along channel. L=20 nm, 
t„=5 nm, Vg = 0.1 V, Vds = O.OIV, VKB.f = VKB.b = 0.0 V 



Figure 2.12. Effect of silicon body thickness variation on surface potential along channel. 

L=50nm, t„,=2 nm, Vg = 0.1 V, Vos = O.OIV, V|rB,f = VirB,b = 0.0 V 
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Chapter 3 


Modeling of Sub-Threshold Current 


In this chapter first an analytical model for sub-threshold current for a long channel 
DG MSOFET will be derived and then it is shown how it can be converted to a model 
for short-channel devices using the expression for the surface potential derived in 
chapter 2. 

Ideally it is expected that when the gate voltage is below the threshold voltage 
there should be no current but in reality the situation is different. Even when gate 
voltage is below the threshold voltage there are some charge carriers, which are very 
small but not exactly zero. So this leads to the flow of some finite drain current even 
when gate voltage is less than threshold voltage. This is what is know as the sub- 
threshold leakage and this constitutes the sub-threshold current. 

It has already been shown with the help of simulation and also on the basis of 
physics based models that diffusion of carriers is the main mechanism of charge 
transport in the sub-threshold region. In sub-threshold region some of the minority 
carriers are able to overcome the barrier at source-channel region and are then swept 
to the drain due to the lateral field and this constitutes the sub-threshold current. 

3.1 Diffusion current 

In a semiconductor at any time the current density at any instant is given by 

Jn ~ n, drift '^n,dijf 
J p = '^p, drift p.diff 

Thus both drift current and diffusion current make up the total current in a 
semiconductor. They may not be occurring at the same time, but the equation is still 
valid. In the sub-threshold region the dominant mechanism is diffusion current as 
mentioned earlier. For a long-chatmel device it has been assumed that the potential 
along the channel is constant which results in a zero electric field thus making 
contribution due to drift current, to sub-threshold current, zero. But as the channel 
length is reduced this assumption of constant surface potential along the channel is no 
longer valid and thus the minimum potential along the channel needs to be found out 
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which will be used for calculating the diffusion current, as at the point of minimum 
potential the electric field will be zero. This results in the total current due to diffusion 
current. Also for a short channel device the point of minimum potential is unique 
unlike long-channel devices where there is a large part along the length where 
potential is constant. Thus in sub-threshold region re-writing the equations for current 
densities 

■^n='^n,diff (3.3) 

Jp = Jp,diff (3-4) 

As an n-DG MOSFET will be considered so the current due to electrons will be 
calculated. Similar formulations can be obtained for a p-DG MOSFET by replacing n 
with p and q with -q. Diffusion current can be expressed as 

dn 

J n,diff - (3-5) 

where Z)„ = diffusion constant and, 

— = concentration gradient of minority carriers 
dy 


3.2 Long Channel Current Formulation 


From Eq. (3.5) the expression for sub-threshold current is given as [17] 

(3.6) 

dy 

In present case there is volume inversion, as reported in [20], hence the complete 
silicon body thickness is taken. Thus the current flows through whole of the silicon 
body. 


If Vchanneliy) is the potential variation along channel direction due to applied 
drain bias, then the expression for the carrier concentration variation with distance 
from source is given as 


n{y) = rii exp 


Cy ) ^channel (y ) 

<t>t 


(3.7) 


Now for a long-channel device the surface potential along the channel is assumed to 
be constant with value 
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(3.8) 


Vs - ~ ^FB 

Thus this results in modification of Eq. (3.7) to 
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where — = Thermal Voltage 
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Thus using Eq. (3.9) in Eq. (3.6) and integrating from source to drain 
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(3.10) 


Also Einstein’s relation relates diffusion coefficient to the mobility as follows 
Dn=<f>tl^n (3.11) 

Thus expression for a long chaimel device, sub-threshold current is 

(3.12) 


W 


^DS = 


7 

1 ■ 

1 - exp 


V j 


^ )_ 


3.3 Short Channel Current Formulation 


Eq. (3.12) will work fine for the devices where surface potential is constant 
for a large part of channel length. But as the length is decreased further, the potential 
does not remain constant along the channel, which is the underlying assumption made 
in the derivation of Eq. (3.12). Since potential is not constant so minimum potential 
along the channel is found at the point where electric field will be zero. This will 
result in zero contribution due to drift current. For evaluating minimum surface 
potential, Eq. (2.9) is differentiated and equated to zero to give the point along the 
channel where the surface potential will be minimum. The minimum potential point 
along the channel (ymm) is given by 
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(3.13) 


Now using this value of and evaluating surface potential at this point using 
Eq. (2.9), the value of minimum surface potential along the channel is given by 

%min = ^Gejf + ^^^|^<^2i<Ps(o) " VoeffhsiL) “ (3.14) 


where 
Cj = sinhl 

Cl = cosh 


'L 

X 

(D 


and 


It can be seen by calculations that for a long-channel device where (L » X) Eq. (3.14) 
reduces to Eq. (3.8) thus resulting in a unification of surface potential for long channel 
and short channel device. The derivation for Eq. (3.13) and Eq. (3.14) is given in 
Appendix C. 

Diffusion current is evaluated at minimum potential point along channel, as is 
done by Liang et.al [1 1], as sub-threshold conduction is determined by the minimum 
potential along the channel which would produce the maximum barrier to carrier 
flow. Thus to get expression for short-channel sub-threshold current the value of 
surface potential of long channel is replaced with the minimum surface potential 
expression as given in Eq. (3.14). Now looking at Eq. (3.12), the surface potential for 
a long-channel is given by {Vc-V^b) ^enn. Thus replacing this with <Ps min 
following expression for sub-threshold current for a short-channel device is obtained 


a ^ I 

^DS ~ 


f \ 

min 

1 - exp 

(-Vds)^ 

1 J 


^ J. 


(3.15) 


3.3 Sub-threshold Swing Formulation 

One of the important performance measurement parameter for a MOSFET in 
sub-threshold region is amount of gate voltage swing required to switch device from 
its ON state to OFF state. This gate voltage swing is called inverse Sub-threshold 
slope (S). It is defined as the amount of change in gate voltage required for producing 



one-decade change in sub-threshold current. So for achieving better switching 
characteristics it should be small as possible. Clearly, S is a convenient measure for 
turn-off characteristics of a MOSFET. By definition it is given as [ 17 ] 




= 2.303 




{mV I decade) 


(3.16) 


Now from Eq. (3.15) 

<i'(ln I ^ d(pg ^ 1 d.(p^ 
dVg dVg (j)^ dVQgjjf 


Using Eq. (3.14) one gets 

min ^ J ^ ^^2 k (Q) ~ ^Geff + <Ps (^) " ^Geff J 

^Ci^^C2(<Ps{0)-VGeffl<Ps{L)-VGeff)- 


C2[Ps{^)^(Ps{L)-2VGeff\ 

Cl ^2C2 (<Ps (0) - VGeff h (h - Voeff ) “ vk 


(3.18) 


Thus final sub-threshold swing expression for symmetric DG MOSFET is given as 

S= 2303<p, 1 + 

In Eq. (3.19) it can be seen that the sub-threshold swing is dependent upon 
through . So to make sure that device is well within sub-threshold region and for 
comparison with the results obtained with 2-dimensional device simulation is kept 
as 10^, 


C2[<Psi0)+<Ps{i^)-2VGejr] 
ClJ2C2{<Ps (O) - VGejf \9s i^) “ ) “ ^DS 


(3.19) 


3.5 Results of comparison with 2-D numerical simulation 

The basic device structure shown in chapter 2 is made in 2-D device simulator 
and simulated using Atlas, which solves the basic device equations like Poisson’s 
equation, continuity equation at the grid points created in simulator by discretizing 
these equations and solving iteratively for various bias points at each of the grid point. 
Mid-gap gates^ are considered in both numerical simulation and analytical model 
calculations. Results have been verified by varying channel length from 300 nm down 


^ Refer Appendix A 
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to 10 nm, silicon thickness from 5 nm to 20 nm in steps of 5 nm and 4 oxide thickness 
values ot 1.0 nm, 1.2 nm, 1.5 nm and 2.0 nm are used. Along with that constant 
mobility model is used with mobility value of 1400 cm^V^sec’^ 

3.5 J Current Comparison Results 

Figures 3,1 - 3.4 show the comparison of model results with numerical simulations 
for different oxide thickness and silicon body thickness with a maximum error of 20 % 
in case of short-channel devices. The normalized sub-threshold current (Ids x L/W) 
obtained using analytical model {Eq. 3,15) is compared with 2-D numerical simulator 
[ 12 ] results for various lengths, silicon body thickness and oxide thickness and 
parameter rj values are obtained. Using these values of 7 and exponential regression 
analysis the following expression for 7 is obtained with coefficient of determination 
(COD) value of 0.95. 

^ Jr364) (0.2617) (0.67358) r(~0.88869) 

• ^ox ^ (3.20) 

where L, tsi, tox are in nm. 



Figure 3.1. Comparison of analytical current model Eq. (3.15) with 2-D atlas simulation for 
varying lengths. Vds= 1.0 V, t,i = 5 nm, t„, = 1.5 nm, L = 20, 30 and 50 nm, VFB,f = VFB,b ~ 0 V , 
mobility (p) = 1400 cm*V‘sec'‘ 
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Figure 3,2. Comparison of analytical current model Eq. (3.15) with 2-D atlas simulation for 
varying lengths. Vos“l«0 V, t,i = 10 nm, to* = 1*2 nm, L = 20, 30, 50 and 70 nm, VFB,f “ = 0 V , 

mobility (p) = 1400 cm^V'^sec"^ 



Figure 3.3. Comparison of analytical current model Eq. (3.15) with 2-D atlas simulation for 

• t .. */ _■! A v . — * 1 ^ nm t 1.0 nm. L “ 30, 50. 70 and 100 nm. 

varying lengths. Vds= 1.0 V, t., - 13 nm, i., i.w > ’ 

VpB,f « VpB b “ 0 V , mobility (p) = 1400 cm Vsec ' 
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Figure 3,4. Comparison of analytical current model Eq. (3.15) with 2-D atlas simulation for 
varying lengths, Vds* 1*0 V, t,i = 20 nm, to* - 2,0 nm, L = 50, 70, 150 and 300 nm, 
VFBvf ** ^FB,b “ 0 V , mobility (|i) * 1400 cm^V’*sec'* 

3.5,2 Sub-threshold Slope Comparison Results 



Channel Length, L (nm) 

Figure 3.5. Comparison of analytical sub-threshold slope model Eq. (3.19) with 2-D atlas 
simulation for varying lengths. Vds= 1.0 V, t., = 10 nm, C. = 1.5 nm, L = 20, 30, 50, 70, 100, 150, 

200, 250 and 300 nm, VFB.f = VFB.b = 0 V , mobility (p)= 1400 cm^V 'sec 
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Figure 3.6. Comparison of analytical sub-threshold slope model £q. (3.19) with 2-D atlas 
simulation for varying lengths. Vds= 1.0 V, t,i = 15 nm, t,, = 2.0 nm, L = 30, 50, 70, 100, 150, 200, 
250 and 300 nm, VFB,f = Vkbj, == 0 V , mobility (p) = 1400 cm^V’sec'* 

Figure 3.5 and figure 3.6 show the comparison of sub-threshold slope values obtained 
using Eq. (3.19) and those extracted from 2-D Atlas simulation. A good agreement is 
found between the two showing the validity of Eq. (3.19). 

3.5.3 Threshold Voltage Roll-off Comparison 

For comparison of threshold voltage roll-off, the most widely used definition 
of threshold voltage based on constant current is used. As is mentioned in Arora [17], 
current can be chosen between 10'^- 10 ^ A after normalizing it with W/L. In present 
work, the current level is taken as 10'® A. Thus threshold voltage is calculated for 
different lengths at the points where Ids ” W/L* 10® A. For calculating threshold 
voltage roll-off, the longest channel threshold voltage is considered as the reference 
and the decrease in threshold voltage for smaller lengths is calculated with respect to 
this long channel threshold voltage. A good agreement is found between the roll-offs 
obtained using analytical model results and 2-D simulations results. Along with this 
the effect of drain voltage on roll-off is also studied. The results of comparisons are 
shown in figures 3.7- 3.10 below. 
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Figure 3*7. Comparison of threshold voltage roll-off obtained from analytical model Eq. (3.15) 
with 2-D atlas simulation for varying silicon body thichess. V V, to* = 1*5 nm, L = 8, 10, 

20, 30, 50, 70, 100, 150, 200, 250 and 300 nm, VFB.f = VFB,b = 0 V , mobility (n) = 1400 cm^V'sec ' 



Figure 3.8. Comparison of threshold voltage roll-off obtained from analytical model Eq. (3.15) 
with 2-D atlas simulation for varying silicon body thichess. V ds= 0-01 V, t„ = 2.0 nm, L = 8, 10, 
20, 30, 50, 70, 100, 150, 200, 250 and 300 nm, Vfb,= VFB,b = 0 V , mobility (p) = 1400 cmV'sec ' 
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Channel Length, L (nm) 

Figure 3.9 Comparison of threshold voltage roll-off obtained from analytical model Eq. (3.15) 
with 2-D atlas simulation for varying drain bias, showing effect of DBL. t = 10 nm, t„, 1.5 nm, 

L - 20, 30, 50, 70, 100, 150, 200, 250 and 300 nm, VFB,r = Vpe.b = 0 V , 
mobility (p) “ 1400 cmV'sec*' 



Figure 3.10. Comparison of threshold voltage roll-off obtained from analytical model Eq. (3.15) 
with 2-D atlas simulation for varying drain bias, showing effect of DBL. t ,i = 20 nm, t„, = 2.0 nm, 

L = 50,70, 100, 150, 200, 250 and 300 nm,VF8j= VFBj. = 0 V, mobility (p) = 1400 cm^V sec 
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Chapter 4 

Modeling of Surface-Potential for Asymmetric DG MOSFET along the 
channel in the Sub-Threshold Region using Gauss’s Law 


4.1 Gauss Law 



The asymmetric device that is considered in this work is operated in 
asymmetric mode by varying the work function of the top and bottom gates. The 
mode of operation of the device is as illustrated in figure 4.1 below with top gate as 

n+ poly-silicon {V^j = -0.56K) and bottom gate as p+ poly-silicon = 0.56 V ) . 


Figure 4.1. Asymmetric DG MOSFET Structure 

In the device, the same gate voltage is applied to both the gates and the same oxide 
thickness is kept for top and bottom gate. Applying the Gauss’s law to the silicon 
body of the asymmetric DG MOSFET 

'-f (4-1) 

-in ^ 

where the firet teim on the L.H.S represents the net electric flux entering the silicon 
body from the lateral direction and the second and the third term represents the net 


electnc flux entering the body from the top and bottom. The R.H.S represents the net 
charge enclosed by the Gaussian box. 


In case of sub-threshold region, the mobile charge carriers will be neglected. 
Thus Eq. (4.1) becomes 





S<P{x,y) 

dy 


dx-e^^ 


y V 

0 0 


(4.2) 


4.2 Assumption for potential variation in perpendicular 
direction 

It has been shown in the earlier works of Yuan Taur [16] and Gen pei et. al [7] 
and also verified by 2-D Atlas simulations that in the sub-threshold region, potential 
variation along the perpendicular direction (i.e. x-direction) is linear. One of the plots 
showing the variation of potential along x-direction is shown in figure 4.2 below. 



Figure 4.2. Variation of potential perpendicular to channel, t,, = 20 nm, t„i = 2.0 nm, Vg - 0.1 V, 
Vds = 1.5 V, L=300 nm. (Simulation done using Atlas) 

Thus one can assume potential to be varying linearly along the perpendicular direction 
and can be written as 
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^(x,jk) = k^x + k^ 

with the following boundary conditions 


(4.3) 


<P -- =<P,j{y) = Front-Gate surface potential 

V 2 y 


(p (y) = Back-Gate surface potential 

V 2 y 


Based on Eq. (4.4) and Eq. (4.5) it can be shown that ^(x,y) is given by 

, ^ -<»../) ,(<Ps^+<Psj) 

^(x,y) = ; x + 4 ^ 


Differentiating Eq. (4.6) with respect to y one gets 

d<p(x,y} ^xf d(p,^,{y) ^ 1 f ^ 

dy t,,\ dy dy ) 2[ dy dy ) 

Integrating Eq. (4.7) from to 

'“f 1 [ '"f xdx + -( '"f 

Jn ^si[ ^ ^ J-i.n A ^ ^ J-J/2 


\[ ^9Ay) 


^si I ^ 


3;^ P-U2 H ^ 


tJ ^9, Ay) , ^9s.f{y) 


2\ dy dy ) 

As nientioned earlier that for short-channel device, potential no longer varies linearly 
in the x-direction» Thus parameter is introduced to adjust the lateral flux for short 

channel devices, Eq. (4.8) is then modified to 

( ^9sAy) I (4.9) 

-Jn ^ 2;;^ ^ J 


Eq. (4.2) and Eq. (4.9) give 

Ml r ^ 1 + f „ V / y Vy ' + f (y ' = 0 

2 tj dy dy j 5 ' o 


(4.10) 


^ Refer section 5.4. 1 
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Now ^(.v') and £„, ,,{)>') are the fields entering from the top and bottom gates and 
are given by 




\yG-yfBj-<psAy)\ 


(4.11) 




\VG-yFB.B-<P..t^y)\ 


(4.12) 


Putting 


tag Eq. (4. 11) and Eq. (4.12) in Eq. (4.10) and differentiating with respect to y 


Sik 

Irj 


d^<p.,iy) . d'^VsAy) 


d/ 


‘ + 


df 


+ 


f ^ ‘'ox 


Now, by solving 1-D Poisson’s equation one can get, [13] 
O’s.i =(»>■./ - £>■•/ ‘si -fsj -^(>'0 - ''FB.f - ‘fsj) 


<P..r 


'' C 

l+~ 

V '^si J 




"'si 


From Eq. (4.14) one gets 

d^cps.b _ “ I ^£ox 

dy^ dy^ V 2 

Using Eq. (4.14) and Eq. (4.15) in Eq. (4.13) one obtains 

dy^ 

where A 




^si^si^ox. s Characteristic Length , 


^GeffJ 

f 


Vq -aVfBj yFB,b . 


a = 


1 + 


C ^ 


•'SI 


2 + 


■'OX 

^si J 


= Gate coupling factor 


(4.13) 


(4.14) 


(4.15) 


(4.16) 
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% = (0‘x)(l) 


4.3 F ormulation of boundary Conditions 

Now to formulate the boundary conditions an undoped body DG MOSFET is 
assumed with n+ source and drain. 


xa 




Figure 43. Boundary conditions for potential at various surfaces of the silicon body for 

asymmetric DG MOSFET 

The boundary conditions, as shown in figure 4.3, are 

, / (0) = <» (- -y .0) = ® 

<Ps,f (^) = «9(- Y , i) = hi + ^DS 

<Ps,b W = > ■^) = hi + 

4.4 Analytical Model for Surface Potential in Sub- 
threshold Region 

Eq. (4.16) is same as Eq. (2.8) with similar boundary condition. It is solved using 
boundary conditions (I) and (II) as is done for Eq. (2.8) in appendix B to give the 
following expression for the front gate surface potential along the channel in the sub- 

threshold region ■ 





9s, f ^y) ^Geff,f + ^s, f ( 0 ) ~ ^Geff,f ) 


sinh 
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sinh 
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sinh 


sinh 
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(4.17) 


Similar expression for back gate surface potential along the channel can be found 
given in Eq. (4.18) below 


as 


9s,b (j) = ^Gejr,b + (p^.aCO) - VGejf,b) 
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sinh 

(L-y'] 
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^ X J 
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+ [9s,biL)-VQgffi,) 
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sinh 
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sinh 
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> 


where = Vq - a V^s,b " 0 - 


(4.18) 


From Eq. (4.17) and Eq. (4.18) it is clear that both front gate surface potential and 
back gate surface potential show similar behaviour as shown by surface potential 
along channel for symmetric DG MOSFET case. Figure 4.4 and figure 4.5 show the 
calculated results based upon Eq. (4.17) and Eq. (4.18) along with its comparison with 
the 2-D numerical simulations obtained using Atlas device simulator for small drain 
bias. 



Figure 4.4. Comparison of Surface Potential along channel for analytical model Eq. (4.17) and 
Eq. (4.18) with 2-D Numerical Simulation (longer lengths) for small drain bias (Vds = 0.01 V), 
t,, *2.0 nm, t,i * lOnm, Vg *-0.1 V, VFB,f= -0.56 V,VFB,b = 0.56 V, L=100, 150 and 300 nm. 







Figure 4.5. Comparison of Surface Potential along channel for analytical model Eq. (4.17) and 
Eq. (4.18) with 2-D Numerical Simulation (short lengths) for small drain bias (Vos = 0.01 V), 
- 2.0 nm, t„ - 10 nm, Vc ■= -0.1 V, Vre,r = -0.56 V, VpB.b = 0.56 V, L=30, 50 and 70 nm 



Figure 4.6. Comparison of Surface Potential along channel for analytical model Eq. (4.17) and 
Eq. (4.18) with 2-D Numerical Simulation (longer lengths) for large drain bias (VDS = 1.5 V), 
tox «2.0 nm, tsi - 10 nm, VG = .0.1 V, VFB,f =-0.56 V, VFB,b = 0.56 V, L-70, 150 and 300 nm. 








Figure 4.7. Comparison of Surface Potential along channel for analytical model Eq. (4.17) and 
Eq. (4.18) with 2-D Numerical Simulation (short lengths) for large drain bias (Vos = 1-5 V), 
t., = 2.0 nm, t„ - 10 nm, Vc = -0.1 V, VFB,f = -0.56 V, Vfb.s = 0.56 V, L=30, 50 and 70 nm. 

From figure 4.4 and figure 4.5 it is clear that potential shows behaviour similar 
to symmetric DG MSOFET. As the channel length is decreased, potential no longer 
remains constant along the channel. Figure 4.6 and figure 4.7 show the variation of 
surface potential along the channel for larger values of drain voltage indicating the 
validity of the model for higher values of also. 




Chapter 5 


Modeling of Sub-Threshold Current 


In chapter 3 it was discussed that in sub-threshold region the current is mainly due to 
diffusion so only diffusion component of current will be considered for formulation of 
the model. 

5.1 Diffusion current 


Sub-threshold current densities can be expressed as 


dn = dn,diff 

(5.1) 

II 

(5.2) 


Since an n-DG MOSFET is considered so the current due to electrons will be 
calculated. Diffusion current is given by 


r n 

Jn,diff = 


(5.3) 


From Eq. (5.3) the expression for sub-threshold current is given by 


{ tJ2 
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(5.4) 
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where is the potential variation along channel direction due to applied drain 

bias and is a function of y only. 

Mobile charge expression is given by 

<p{x, y) ~ ^channel iy') 


n{x, y) = exp 




Thus one gets, 
dn{^x,y) 


dV.., 


n. exp 


^ channel 

Mng Eq. (5.6) in Eq. (5.4) one gets, 
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channel 
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(5.7) 



Evaluation of «(>;) is done by integrating Eq. (5.5) from the top surface to the 


bottom of silicon body (Appendix D) 

expf^l-expf^l 

n (y) = ■' s . ^ i -^ expf = _^expf 

\(<PsA {•!>,) 3l(y) 1, <j>, 


<t>, ) \<f>, 


where 


iR(y) = 


(Psj, [ (PsJ 

<!>, ) I 


f Kb { <Ps.f 

\<i>, ] I 


Thus using Eq. (5.8) in Eq. (5.7) one gets 

I - aP ~0 expf 1 ^^cAom 

I <!>, J dy 


Integrating from source to drain gives 


Ias=qrip„Wt, 1-exp -f- 

\ <P> 


J9t(y)dy 


(5.10) 


Also Einstein’s relation relates diffusion coefficient to the mobility as follows 
Dn==<t>tMn (5.11) 

Thus expression for sub-threshold current becomes 

= 1-exp j \jfi{y)dy (5.12) 

From Eq. (5.12) it is clear that it can not be solved analytically. So we need to come 
up with some constant K such that the integral in the denominator of Eq. (5.12) can be 


written as 


\'i{{y)dy = K\dy = K.L 


(5.13) 


K is nothing but the average value of integral over the range of integration. Thus, if 
one can come up with some analytical expression for K then Eq. (5.12) can be 
converted into an analytical expression. ■ 
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5.2 Evaluation of Integral K 


Function is plotted in figure 5.11. The average value (K) of 3l(>’) is 

calculated numerically. Also plotted in the figure is the back surface potential as a 
function of y. 
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Figure 5.1. Comparison of 91 (y) with back surface potential along the channel. t.i = 10 nm, 
U = 1.5 nm, Vg * -0.1 V, Vos = 1.5 V, L = 30 nm, Vra,f = -0.56 V, VFB.b = 0.56 V. 

From the figure it is clear that the average value of 91 (y) is occurring close to the 
point of minimum back surface potential (ymin.^;)- So if one evaluates the function 

under integral at the point where the minimum of back surface potential is occurring 
then a reasonable value of K can be obtained. Based on this assumption, K can be 
represented as 


K = 



‘Ps,b ) 

A j 

- 

f 

K A 


exp 

<PsA 

) 

exp 

<l>t JJ 


(5.14) 


y ymin.ft 
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where is the point along channel where the back surface potential is 

minimum. Now be found using the similar procedure as applied for 

symmetric case for finding out the minimum surface potential point. From the 
.procedure outlined in appendix C one gets the following expression for ^ 




^In 


(0) - VGeff,b ) expj^l j - (L) - Voeff^b ) 

^s.b i^) - ^Gejf,b ) - ^s,b (0) - ^Geff,b ) exp 


-L 

A 


(5.15) 


From Eq. (5.12), Eq. (5.14) and Eq. (5.15), the final analytical expression for sub- 
threshold current can be written as 


r .4 ^ 

^ DS ~ ^ ^ ^si 


1 - exp 


' Vds '^ 


(5.16) 


5.3 Sub-threshold Swing Formulation 


As mentioned in chapter 3, sub-threshold swing is an important performance 
measurement parameter for a MOSFET in sub-threshold region. Sub-threshold swing 
is given by 
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(5.17) 


Now from Eq. (5.16) 
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Using Eq. (5.14) one gets 
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(5.18) 


(5.19) 


The mathematical analysis done to obtain Eq. (5.19) is shown in appendix E. Thus 
final sub-threshold slope expression for asymmetric DG MOSFET is given as 
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5 = 2.303 1- 

In Eq. (5.20) it can be seen that the sub-threshold swing is dependent upon Vg 
through y . So to make sure that device is well within sub-threshold region and for 

comparison with the results obtained with 2-dimensional device simulation, Vq is 
kept as -5^1 . 

5.4 Results of comparison with 2-D numerical simulation 

The basic device stmcture shown in chapter 4 is made in 2-D device simulator 
and simulated using Atlas. Results have been verified by varying channel length from 
300 nm down to 20 nm, silicon thickness from 5 nm to 20 nm in steps of 5 nm and for 
4 oxide thickness values of 1.0 nm, 1.2 nm, 1.5 nm and 2.0 nm. Along with that 
constant mobility model is used with mobility value of 1400 cm^V'sec"', 

From figures 4.4 - 4.7 one can see that surface potential is not exactly 
matching with the 2-D simulation results. But as far as predicting the value of K is 
concerned, the proposed model is able to predict it within reasonable error range 
compared to 2-D simulation results as shown in figure 5.2. The squares in the figure 
show the enror as a function of y where fe is given by 

xlOO (5.21) 

Also is shown the variation of back surface potential along the channel. The back 
surface potential is calculated using the analytical expression using Eq. (4.18). 

From the figure it is clear that the percentage error in evaluation of value of 'R (y) 

is minimum in the region close to j as is outlined by the box in figure 5.2. Now K 
is evaluated at ^ . This explains why the approximation gives a reasonable match 
of model with the 2-D simulation results as shown in following sub-sections. 


nyt. 
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Figure 5.2. Variation of percentage error f, along the channel. V^s = 0.01 V, t„i = 2.0 nm, 
t„ - 10 nm, Vc “ -0.1 V, -0.56 V, VFB.b = 0.56 V, L=30 nm. 

5.4.1 Current Comparison Results 

Figures 5.3 — 5.6 show the comparison of model results with numerical 
simulations for different oxide thickness and silicon body thickness. Now due to 
simplicity of model and the simplifying assumptions made in formulations of the 
model, the maximum error in predicting current values have been under 50%. Similar 
to what has been done for symmetric case, the expression of parameter rj found 
using exponential regression analysis with COD value of 0.952, is given as 

( 2 . 01028 ) , ( 0 . 33548 ). ( 0 . 9958 !) ,(- 1 . 07241 ) (5.22) 

” ~ ‘ot 

where tji, tox and L are in nm. 


40 





Figure 5.3. Comparison of analytical current model Eq. (5.16) with 2-D atlas simulation for 
varying lengths. Vds= 1.0 V, t,i = 5 nm, to, = 1.5 nm, L = 20, 30 and 50 nm, = -0.56 V, 
VFB,b - 0-56 V, mobility (p) = 1400 cm^V'sec"' 



Figure 5.4. Comparison of analytical current model Eq. (5.16) with 2-D atlas simulation for 
varying lengths. Vds= 1.0 V, t„ = 10 nm, t,, = 1.2 nm, L = 20, 30, 50 and 70 nm, VFB.f = -0.56 V, 
Vf.B.b = 0.56 V, mobility (p) = 1400 cm Vsec*' 
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Figure 5.5. Comparison of analytical current model Eq. (5.16) with 2-D atlas simulation for 
varying lengths. Vds=1-0 V, t,i = 15 nm, to* = 1.0 nm, L == 30, 50, 70 and 100 nm, VFB,f = -0.56 V, 
VpB^b = 0.56 V, mobility (p) = 1400 cmV"*sec'* 



Figure 5.6. Comparison of analytical current model Eq. (5.16) with 2-D atlas simulation for 
varying lengths. Vos— 1*0 V, tsj = 20 nm, to* = 2.0 nm, L = 50, 70, 150 and 300 nm, VpB,f ~ -0.56 V, 
VFB,b = 0.56 V, mobility (p) = 1400 cm^V 'sec"' 




5.4.2 Sub-threshold Swing Comparison Results 



Figure 5.7. Comparison of analytical sub-threshold swing model Eq. (5.20) with 2-D atlas 
simulation for varying lengths. Vds=1.0 V, t,i = 10 nm, to, = 1.5 nm, L = 20, 30, 50, 70, 100, 150, 
200, 250 and 300 nm, VFB,f = -0.56 V, Vre.b = 0.56 V, mobility (p) = 1400 cmVsec ' 



Figure 5.8. Comparison of analytical sub-threshold swing model Eq. (5.20) with 2-D atlas 
simulation for varying lengths. Vds= 1'0 V, t,i = 15 nm, t„i = 1.2 nm, L = 30, 50, 70, 100, 150, 200, 
250 and 300 nm, VpB.r = -0.56 V, VpB,b = 0.56 V, mobility (p) = 1400 cm V sec ' 


Figure 5.7 and figure 5.8 show the comparison of sub-threshold slope values obtained 
using Eq. (5.20) and those extracted from 2-D Atlas simulation. A good agreement is 
found between the two showing the validity of Eq. (5.20). 

5.4.3 Threshold Voltage Roll-off Comparison 

As is done for symmetric DG MOSFET, comparison of threshold voltage roll-off is 
done for asymmetric DG MOSFET. A good agreement is found between the roll-offs 
obtained using analytical model results and 2-D simulations results. Along with this 
the effect of drain voltage on roll-off is also studied. The results of comparisons are 
shown in figures 5.9-5.12 below. Here also the current is normalized with respect to 
W/L and current level is taken as 10'® A after normalization for defining the threshold 
voltage of the device. 



Channel Length, L (nm) 


Figure 5.9. Comparison of threshold voltage roll-off obtained from analytical model Eq. (5.16) 
\vith 2-D atlas simulation for varying silicon body thickness. Vds*1»0 V, to, =* 2.0 nm, L * 20, 30, 
50, 70, 100, 150, 200, 250 and 300 nm, VpB^f == -0.56 V, == 0.56 V, 

mobility (p) = 1400 cm^V'*sec"* 
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Figure 5.10. Comparison of threshold voltage roll-off obtained from analytical model Eq. (5.16) 
with 2-D atlas simulation for varying silicon body thickness. Vds= 0.01 V, to* = 1.5 nm, L = 20, 30, 
50, 70, 100, 150, 200, 250 and 300 nm, VpB.r = -0.56 V, VpB.b = 0.56 V, 

mobility (p) = 1400 cm^V"*sec'‘ 
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Figure 5.11. Comparison of threshold voltage roll-off obtained from analytical model Eq. (5,16) 
with 2-D atlas simulation for varying drain bias, showing effect of DIBL. Lj = 10 nm, tox = 1.5 nm, 
L = 20, 30, 50, 70, 100, 150, 200, 250 and 300 nm, VpB^f = -0.56 V, VFB.b = 0.56 V, 
mobility (p) = 1400 cmV'sec ' 
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Figure 5.12. Comparison of threshold voltage roll-off obtained from analytical model Eq. (5.16) 
with 2-D atlas simulation for varying drain bias, showing effect of DIBL. t,i = 20 nm, t„, = 2.0 nm, 
L “ 50, 70, 100, 150, 200, 250 and 300 nm, VpB.r = -0.56 V, Vpg.b = 0.56 V, 
mobility (p) » 1400 cm^V'sec’' 
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Chapter 6 

Symmetric Vs Asymmetric DG MOSFET - A Comparative analysis based 
on the proposed model 


In the previous chapters, models for sub-threshold region for asymmetric and 
symmetric DG MOSFETs have been presented. There comparison has been done with 
2-D device simulator and a good agreement between the two is found showing the 
validity of models made. But the whole analysis is not complete until a comparison is 
done between the two modes of operation of these devices and figure out how one 
mode is better than the other mode. A lot of wok has been done in comparison of 
these two modes of operation [17],[20]. Most of the work point towards the 
usefulness of asymmetric devices towards threshold voltage control keeping the 
similar roll-off values. In this chapter a comparative analysis will be done for both 
modes based upon the comparison of sub-threshold current, threshold voltage roll-off, 
the effect of drain induced barrier lowering and sub-threshold slope by comparing two 
devices of same parameters except change in the gate work function. Finally an 
attempt is made to come up with a unified model for DG MOSFETs, applicable to 
both symmetric and asymmetric DG MOSFETs 

6.1 Sub-threshold Current Comparison 

Figures 6.1 and 6.2 show the comparison of sub-threshold currents of 
symmetric and asymmetric DG MOSFETs with only difference that for symmetric 
DG MOSFET, mid-gap gates are considered and for asymmetric DG MOSFET, n+ 
and p+ gates are considered and all other device parameters are considered same. 
From the figures it is clear that with same device parameters asymmetric DG 
MOSFET has lower threshold voltage compared to symmetric DG MOSFET as is 
also expected as in asymmetric case one gate is n-t- so this will have lower threshold 
voltage and will turn on early. Similar conclusions have been drawn by Yuan Taur 
[16], where he has commented on the tum-on behaviour of symmetric and 
asymmetric DG MOSFETs on the basis of inversion charge under the gate in the two 
cases considered. Thus the results obtained are consistent to what has already been 
shown by analyzing it in terms of current. 
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Gate Voltage, (V) 


Figure 6.1 : Comparison of symmetric and asymmetric, analytical current model with 2-D atlas 
simulation. Vds=1-0 V, t,i = 20 nm, t^x = 2.0 nm, L = 50 nm, Vfbj = VFB,b == 0 V (symmetric), 
VFB,f ==-0.56 V, VFB,b - 0.56 V (asymmetric), mobility (p) = 1400 cm^V^sec * 




Gate Voltage, V^ (V) 


Figure 6.2 : Comparison of symmetric and asymmetric, analytical current model with 2-D atlas 
simulation. Vi>s=1.0 V, t^j = 5 nm, t©* ~ 1*0 nm, L = 30 nm, VFB,f ~ ^FB,b — 0 V (symmetric), 
VFB,r =-0.56 V, VFB,b = 0.56 V (asymmetric), mobility (n) = 1400 cm^V'sec ' 
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6.2 Threshold Voltage Roll-off 


DG MOSFETs have shown good short-channel immunity due to more gate control. 
Now as we go on decreasing the silicon body thickness, the gate control is expected to 
increase on the silicon body and thus resulting in lesser short-channel effects for same 
lengths. This is what has been observed both for symmetric and asymmetric DG 
MOSFETs that as one goes on reducing the silicon body thickness, the threshold 
voltage roll-off reduces for the same channel length. Rephrasing in different way, it 
can be said that a particular threshold voltage roll-off will occur at larger lengths itself 
for larger silicon body thickness. This effect has been observed both for symmetric 
and asymmetric in chapter-3 and chapter-5 respectively. A comparison of roll-off is 
done between symmetric and asymmetric to see which is more immune to short- 
channel effects as shown in figures 6.3 and 6.4 below. 



Channel Length, L (nm) 


Figure 6.3 : Comparison of symmetric and asymmetric, threshold voltage roll-off. V 0 s=l.O V, 
tsi = 5 nm, t„x = 1.5 nm, VpB.r = VFB,b = 0 V (symmetric), VpB.r =-0.56 V, Vps.b = 0.56 V 
(asymmetric), mobility (p) = 1400 cm^V^'sec"'. Results obtained using the proposed analytical 


model. 





Figure 6,4 ; Comparison of symmetric and asymmetric, threshold voltage roll-off. Vds-1-0 V, 
t,i = 15 nm, tox = 2.0 nm, VpB^f = VpB^b = 0 V (symmetric), Vfbj =-0.56 V, VpB.b = 0.56 V 
(asymmetric), mobility (p) = 1400 cmV^sec Results obtained using the proposed analytical 
model. 

Thus from above two figures one can conclude that both symmetric and asymmetric 
DG MOSFETs give similar threshold voltage roll-off but asymmetric threshold 
voltage can be controlled by varying the work functions of the front and back gates. 
This results in better control over device performance. 

6.3 Drain Induced Barrier Lowering 

As we increase the drain bias, this effectively results in the decrease in the 
barrier to the flow of carriers from source to drain and this effect is more prominent 
for short channel devices. This results in lowering of the threshold voltage of the 
device. For comparison of the effect of drain bias on the threshold voltage, the parallel 
shift between the threshold voltage roll-off at Vds = 0.01 V and Vds = 1.0 V is 
calculated for both symmetric and asymmetric devices and is as shown in figures 6.5 

and 6.6. 
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Figure 6.5 : Comparison of symmetric and asymmetric, DIBL effect, VDS=0.01 and 1.0 V, 
tsi = 10 nm, tox = 1.5 nm, VFB,f = VFB,b = 0 V (symmetric), VFB,f =-0.56 V, VFB,b = 0.56 V 
(asymmetric), mpbility (p) = 1400 cm2V-lsec-l. Results obtained using the proposed analytical 
model. 



Figure 6.6 : Comparison of symmetric and asymmetric, DIBL effect, Vds=0,01 and 1.0 V, 
t,i = 20 nm, t„, = 2.0 nm, VpB.f = VFB,b = 0 V (symmetric), VpB.r =-0.56 V, VpB.b = 0.56 V 
(asymmetric), mobility (n) = 1400 cmVsec '. Results obtained using the proposed analytical 
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model. 





Thus from figures 6.5 and 6.6 one can conclude that effect of DIBL is smaller in 
asynunetric compared to symmetric DG MSOFETs. Similar results have been 
reported by J.G. Possum [22], showing the validity of the results obtained. 

6.4 Sub-threshold Swing 

A comparative analysis of the sub-threshold slope obtained for symmetric and 
asymmetric DG MOSFETs will be done in this section. Compared slopes were 
extracted from the 2-D simulation results and also from the sub-threshold swing 
formulation for both symmetric and asymmetric in chapter-3 and chapter-5 
respectively. 



Channel Length, L (nm) 


Figure 6.7 : Comparison of symmetric and asymmetric, sub-threshold slope, Vds= 1-0 V, 
t,i = 10 nm, tox = 1.5 nm, VFB,r = VpB.b = 0 V (symmetric), VpBj —0.56 V, VpB.b = 0-56 V 
(asymmetric), mobility (p) = 1400 cm^V'sec'*. Analytical expression used - Eq. (3.19) 
(symmetric) and Eq. (5.17) (asymmetric). 
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Figure 6.8 : Comparison of symmetric and asymmetric, sub-threshold slope, Vds=1-0 V, 
t„ 15 nm, tox ^ 1.2 nm, VKB,f = VFB,b = 0 V (symmetric), VpB.f =-0.56 V, VpB,b = 0.56 V 
(asymmetric), mobility (p) = 1400 cmV^sec'. Analytical expression used - Eq. (3.19) 
(symmetric) and Eq. (5.17) (asymmetric). 

Thus from figure 6.7 & 6.8 one can conclude that both symmetric and asymmetric DG 
MOSFETs have similar immunity to short channel effects as is indicated by the sub- 
threshold slope variation with channel length. 

All the comparisons done till now show that, in terms of threshold voltage 
roll-off, DIBL and sub-threshold slope both asymmetric mode and symmetric mode of 
operation of DG MOSFETs show similar behaviour. But if one looks carefully, then 
asymmetric DG MOSFET has one main advantage over symmetric and that is the 
dynamic control over threshold voltage. Because of all these things asymmetric DG 
MOSFETs have come up as a viable solution for practical applications and are being 
extensively studied all over the world. 





6.5 Unified Sub-threshold Model 


The similar behaviour obtained for symmetric and asymmetric devices points that 
there is a similarity between the two devices and hence an attempt is made to come up 
with a unified model for DG MOSFETs applicable to both symmetric and 
asymmetric. Looking carefully at Eq. (3.15) and Eq. (5.16) (repeated below in Eq. 
(6.1) and Eq. (6.2) for convenience), then both the equations are same except one 
factor represented by function f{<p,) inEq. (6.3) below. 
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Here K is given by Eq. (5.14) which is 
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Thus one can write a generalized equation for symmetric and asymmetric DG 
MOSFET 
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where f is given by Eq. (6.3) in respective cases. 



Conclusion and Future Work 


Analytical model for sub-threshold region for symmetric and asymmetric DG 
MOSFETs have been proposed. The models have been verified using 2-D numerical 
simulator. The unification of models for symmetric and asymmetric DG MOSFETs 
has been achieved through a surface potential dependent function / ) . The strength 

of the proposed models is that they are physics based, simple and yet are able to 
predict the device behaviour reasonably accurately. They have potential applications 
in circuit simulators. 

Future Work 

DG MOSFETs are in the early stages of development. One of the problems that still 
remains unresolved is a definition of threshold voltage. This device being used with 
un-doped silicon body, hence the classical definition of threshold voltage is not 
applicable to these devices. One has to thus come up with other ways of defining 
threshold voltage especially for un-doped body devices. There are many definitions 
available in literature for defining threshold voltage like constant current, maximum 
slope of trans-conductance etc [23]. But the inherent problem with a definition is 
either have to define threshold voltage at some arbitrary current level or the 
introduction of noise in extracting the threshold voltage by using derivative methods. 
In recent times a new definition for threshold voltage for un-doped body MOSFET 
has been given which defines the threshold voltage as the point of intersection of the 
two asymptotes to linear and saturation regions of operation of device [10]. So there is 
a scope for research in this area. 



Appendix A 

Derivation for comparison of Surface and Centre Potential 


In the derivation of the equations the band structure along with the co-ordinates are as 
defined in figure A-i shown below 



(a) i 


Figure A-1. Schematic band diagram for a symmetric DG MOSFET. (a) At zero gate voltage, 
assuming near mid-gap gates, (b) Band bending when gate voltage is near threshold voltage. 


Since silicon body is not connected to any bias so its Fermi level is referenced with 
respect to the Fermi level of source Fermi level. Hence mid-gap gates means flat band 
voltage is zero with respect to that Fermi level of the silicon body. 

Now applying Poisson’s equation for the silicon body considering only the 
mobile charge carriers as the body is assumed undoped one gets 
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Eq. (A-1) also assumes that q(j)lkT » 1 so one can neglect the charge due to holes 
also. Now multiplying both sides of Eq. (A-1) with 2~ 
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Integrating Eq. (A-2) from any point x in the silicon body till the centre i.e. at x = 0 as 
shown in figure A-1 with the following boundary conditions 
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Rearranging the terms 
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Now integrating Eq. (A-4) again 
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Still simplifying 
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Now surface potential (p^ s <p{xj . Thus from Eq. (A-3) and Eq. (A-6) one gets 
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Now surface potential is related to the gate voltage Vq and oxide thickness tox at Si- 
Si02 interface by Gauss’s law given as 
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Now substituting, 6 ■ 
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and using Eq. (A-3) and Eq. (A-6) 


d(p^ im. 


dx 




exp 


^ me, ure ^ 

2kT 


tan(^) 


(A.IO) 


Also Eq. (A-8) becomes 
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Thus from Eq. (A.9), Eq. (A.IO) and Eq. (A.l 1) one gets 
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Eq. (A.12) and Eq. (A.l I) now form Eq. (2.3) and Eq. (2.4) respectively in chapter 2. 
Eq. (A.12) can be solved using any of the commonly known numerical methods. For 
the present case Newton-Raphson method has been used for solving it. 



Appendix B 

Derivation of Surface Potential expression along the channel for 
Symmetric DG MOSFETs 


From Eq. (2.8) one has 
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Also the following boundary conditions are there from section 2.3 
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Thus Eq. (B.7) forms- Eq. (2.9) in chapter 2. 



Appendix C 

Minimum Surface Potential along the channel derivation for Symmetric 
DG MOSFETs 


Minimum Surface potential point: 
From Eq. (2.9) one has 
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Now differentiating and equating it to zero will give point along channel (ymin) where 
surface potential will be minimum 
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Thus Eq. (C.2) is giving Eq. (3.13) in chapter 3. 
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Minimum surface potential: 


Minimum surface potential along the channel is obtained by evaluating Eq. (CT) at 
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From Eq. (C.3) and Eq. (C.4) one can get 
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Similarly using Eq. (C.5) and Eq. (C.6) one can show 
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Thus Eq. (3.14) in chapter 3 represents Eq. (C.9). 



Appendix D 

Integration of carrier density perpendicular to channel 

Starting with Eq. (5.5) 
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Now integrating Eq. (D.3) from top surface to bottom surface of silicon body. 
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D.4) one gets 
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Appendix E 


Mathematical Analysis for asymmetric DG MOSFET sub-threshold 
slope 


From Eq. (5. 1 8) one has 
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Now from Eq. (4.17) and Eq. (4.18) one gets 


fsjiy) = ^GeffJ + (<^5,/(0) - ^GejfJ ) 


sinh 


L-y 


sinhl ^ 

.2 




sinhf^ 


(E.3) 

li 


sinh, 

V u 




JJ 


(E.4) 


(pa^hiy) = ^GefJ\h + ^s,h (0) - ^Geff,b ) 




I 2 


r j\ 


sinh 


u. 




[9s,bi.i^)-^Geff,b) 


sinh 


sinh 




(E.5) 


69 




Differentiating Eq. (E.4) with respect to Vg, one obtains 



Similarly by differentiating Eq. (E.5), one gets 



From Eq. (E.6) and Eq. (E.7) one can infer 
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Putting Eq. (E.9) in Eq. (E.l) 
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With the help of Eq. (E.8) in Eq. (E.IO) one finally gets 



Thus Eq. (E.l 1) is used as Eq. (5.19) in chapter 5. 
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